Luttinger theorem for a spin-density-wave state 
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We obtained the analog of the Luttinger relation for a commensurate spin-density-wave state. We 
show that while the relation between the area of the occupied states and the density of particles 
gets modified in a simple and predictable way when the system becomes ordered, a perturbative 
consideration of the Luttinger theorem does not work due to the presence of an anomaly similar to 
the chiral anomaly in quantum electrodynamics. 
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One of the key open problems in the study of high- 
T c superconductors is the shape and the area of the 
electronic Fermi surface near half-filling. Photoemission 
studies |lJ performed near optimal doping have demon- 
strated that the area enclosed by the Fermi surface is 
large, and for a doping concentration x, constitutes a 
fraction (1 — x)/2 of the area of the Brillouin zone, just 
as for free fermions. This is consistent with the Luttinger 
theorem which states that in any conventional Fermi liq- 
uid, the area enclosed by the Fermi surface does not 
change due to an interaction between fermions. On the 
other hand, similar experiments on underdoped materials 
have shown no evidence of a Fermi surface crossing for 
momenta far from the zone diagonal Q . These data are 
consistent with the idea that the Fermi surface evolves 
with decreasing doping towards a small hole Fermi sur- 
face consisting of four pockets centered at (±7r/2, ±7r/2). 
The existence of such a Fermi surface in the paramagnetic 
phase would imply a violation of the Luttinger theorem. 

To consider a possible violation of the Luttinger theo- 
rem in a system with strong magnetic fluctuations, one 
first needs to understand how this theorem is modified 
when the system acquires a long-range magnetic order. 
This is the issue which we address in this paper. We 
show that there are several subtleties which emerge al- 
ready on the mean-field level. In particular, a naive per- 
turbative approach yields incorrect results because of a 
hidden anomaly which requires a proper regularization. 

We begin with a brief consideration of the Luttinger 
theorem for a conventional Fermi liquid || . The number 
of particles N can be expressed through the single par- 
ticle Green function Gk{<-*j) of the interacting electrons 
as 
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where Tr implies a summation over spin projections. 
Using the obvious relation for G and the self-energy S, 
Gk{oj) = d log G^ 1 (ui)/duj + Gfc(w)<9£ fe (o;)/<9u;, one can 
rewrite Eq.([l]) as N = Ii — I 2 where 
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In Ref. §, Luttinger and Ward (LW) argued that 
72 = to all orders in perturbation theory around free 
fermions. Their consideration is based on a functional 
Y = -i Tr Y,k I du)/{2ir) Y k (uj) determined by the 
variational equation 
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This variational equation allows one to re-express I2 as 
an integral of a full derivative: 



h = -i Tr 



k 



did dYk(to) 
(27r) du 



(4) 



LW gave a recipe how to obtain the functional Yfc(w) or- 
der by order in a diagrammatic perturbation theory. The 
perturbative Yk(u>) is free from singularities and vanishes 
at large frequencies. Obviously, in this situation I2 — 0. 
(We however will demonstrate below that I2 may be fi- 
nite due to nonperturbative effects). The term I\ also 
contains a full derivative over frequency and, apparently, 
should vanish for similar reasons. However, its calcula- 
tion requires care as log G is a non-analytic function of 
frequency in the upper half-plane. A careful considera- 
tion shows that due to a nonanalyticity, I\ turns out to 
be finite: 
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where [i is the chemical potential of the interacting elec- 
trons. It then follows that 



AT = / 1= Tr ^9(£ fe -/i), 
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where Ek is the quasiparticle energy, and the O-function 
counts the states inside the Fermi surface. This is the 
Luttinger theorem for a conventional Fermi liquid. It 
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states that the number of states enclosed by the Fermi 
surface is equal to the number of particles N. 

The Luttinger relation can be easily extended to the 
case when the system develops a long-range magnetic or- 
der. Consider for defmitcness a 2D spin-fermion model 
which describes a system of propagating fermions with 
dispersion tk coupled to localized spins S q by an exchange 

interaction j§[|: H int = 9E q ,k,a, ( 3 c l,fe+? f ^ c ft' $r 
Here a % are the Pauli matrices, and the localized spins are 
described by their dynamical spin susceptibility, x(g,a;). 
For simplicity, we assume that there are no other in- 
teractions, i.e., for g = one deals with free fermions. 
In the antiferromagnetically ordered phase one of the 
components of S acquires a nonzero expectation value 
< S* >=< S z > 8 q> Q. In the mean-field spin-density- 
wave (SDW) theory, the interaction term reduces to 
^ c a k+Q a aa c a,k, where A = g < S z >. Problems of this 
kind are most conveniently described in a matrix formal- 
ism. In the presence of the SDW order, the electronic 
spectrum can be obtained from the 4x4 matrix 



lj - e k A a z 
A a z us- e k+Q 
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where e k — £fc — /i, and /i is the actual, A-dependent 
chemical potential. For the purpose of comparison 
with the paramagnetic phase we will keep using the 
original Brillouin zone rather than the magnetic one. 
Then, for each k one gets two solutions E k — (Ek + 
e k+Q )/2 ± [A 2 + ((e fe - lk±p)/2) 2 ] 1/2 which satisfy 



E k = E k+ q . The relations (yj),(@) between N and the 
single particle Green's function still hold in the ordered 
state, with the only modification that Tr now acts on 
4x4 matrices. In a mean-field approximation the I2 term 
vanishes identically, and hence N = I±. A substitution 
of Eq.g) into Eq.(§) yields N = J2 k (0(E+) + 0(E k )) . 
Subtracting this relation from the total number of states 
in the Brillouin zone, Nbz, one finally obtains 
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where x = 1 — N/Nbz is a hole doping concentration, 
while Shoie and S double are fractions of the Brillouin zone 
occupied by hole pockets (E ± > 0) and doubly occupied 
electronic states (i? ± < 0), respectively. Eq.(||) is the 
Luttinger relation for the ordered SDW phase. It has 
been derived here in a mean-field approximation, but 
Eq.(||) should also hold in the presence of fluctuations 
due to an interaction between fermions, simply because 
the effect of the spin ordering is completely absorbed into 
the matrix formulation of the LW functional Y. 

We now give the physical interpretation of the Lut- 
tinger relation for the ordered state. For this, we fix the 
doping concentration at some small but finite level and 
follow the Fermi surface evolution with increasing A. For 
A = 0, the Fermi surface has the form shown in Fig |l|a 
- it is centered at (tt, it) and encloses an area which is 
slightly larger than half of the Brillouin zone, in accor- 
dance with the conventional Luttinger theorem, Eq. m). 
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FIG. 1. Schematic evolution of the Fermi surface with the 
gap A in the ordered SDW phase with electron hopping to 
nearest (t) and next-nearest(t') neighbors on a square lattice. 
We set t' = — 0.2t. The dashed line indicates the boundary of 
the magnetic Brillouin zone. 



Switching on an infinitcsimally small A doubles the unit 
cell in real space. As a result, extra pieces of the Fermi 
surface appear because a Fermi surface crossing at k nec- 
essarily implies one at k + Q (Fig. ||b) . The subsequent 
evolution of the Fermi surface proceeds as is shown in 
Fig. |l|c,d. For sufficiently large A, the doubly occupied 
electron states disappear, and the Fermi surface consists 
of just four small hole pockets centered at (tt/2, ir/2) (see 
Fig. |l|d) . The total area enclosed by the pockets is small 
and scales as x in agreement with Eq.(||). 

So far, the LW argumentation works perfectly. How- 
ever, the matrix formalism cannot be extended to the 
paramagnetic phase. We now discuss another, compli- 
mentary approach to the ordered state in which one 
does not explicitly introduce a condensate of S z but in- 
stead obtains the SDW form of the electronic spectrum 
by introducing an exchange of S— functional longitudi- 
nal spin fluctuations: xi(q^) — (A/g) 2 <5(£l)<5(<7 — Q). 
This approach is advantageous as it can be applied to 
the paramagnetic phase - one just has to substitute the 
5— functional longitudinal susceptibility by the one with 
a finite spin-spin correlation length In the spin- 

fluctuation formalism, the full Green's function does not 
have a matrix structure, but there appears a self-energy 
due to spin-fluctuation exchange. To the lowest order in 
A 2 , 



S fc H = A 2 Gi°] c 
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where G k (u>) — (tu — ek)^ 1 is the Green's function of free 
electrons. For the S— functional form of \i this is in fact 
an exact result - the higher order terms in E are absent 
due to a Ward identity [pi. The full Green's function 
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then takes the form 



G fc (w) 



£ k+Q 
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Not surprisingly, it coincides with the diagonal part of 
the matrix G defined by Eq.(Q). At the same time, the 
I\ t 2 terms are not the same in the matrix and the di- 
agonal formalisms. Moreover, the self-energy in (^) has 
a singular frequency dependence. We therefore have to 
reexamine the argument that 1% = 0. To begin with, 
we first observe that I\ alone does not yield the right 
expression for N. Indeed, it follows from Eqs.(]2],|l0|) that 



h= 2J2 (©(4) + 0(4") - 0(^-+q)) 



(11) 



For large A, I\ — 2(1 — x)Nbz in clear disagreement with 
Eq.(||). This implies that I2 should be finite. Observe 
that the numerator in ( |l0| ) also contributes to 1\ . 

To calculate I2, we use Eq.(^) which relates 1% to the 
LW functional Y. This functional can be straightfor- 
wardly obtained by performing the functional integra- 
tion of the self-energy E over G k (ui). To do this, we 
first have to re-express E in terms of the full Green's 
functions. This can be done either directly, by sub- 
stituting Gj^_q(lj) in Eq.(|9|) in terms of G, or dia- 
grammatically, by rewriting the self-energy as T, k ((j) = 
A 2 Gk+Q(oj)^k,k+Q(^>), and collecting diagrams for the 
full vertex T order by order in the perturbative expan- 
sion, each time using full propagators G for the internal 
fermionic lines. It turns out, however, that the results 
of the two approaches are different, because there exists 
a range of frequencies where the perturbative expansion 
is not convergent. Indeed, a formal summation of the 
perturbation series for the self-energy yields E^ er *(w) = 

2A 2 G fc+Q (w)/ (l + v /l + 4A 2 G fc (w)G fc+Q H) . On the 

other hand, the use of an exact relation between bare and 
full Green's functions yields 



2A 2 G fe+Q H 



1± v /l + 4A 2 G fc ( W )G fc+Q H 
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Here the upper and lower signs should be used when 
\y\ < 1 and \y\ > 1, respectively, where y — 

& 2 Gk M^k+Q M • Clearly, the perturbative and exact 
expressions for the self-energy coincide only for \y\ < 1, 
i.e., when the perturbative expansion is convergent. 

We now show that due to the difference between S and 
E pert , the actual I2 is finite, while the perturbative ap- 
proach yields 1% = to all orders in A. Indeed, perform- 
ing the integration over the Green's function in Eq.(^) 
and using Eq.(|l2]) for E, we find 



logG fc H ± ^l + 4A 2 G k {to)G k+Q (uj) 



± l lo ( y/l + 4A 2 G fc ^)G fc+Q n-l 
2 ° g I \/l + 4A 2 G fc (w)G fc+ Q(w) + 1 



(13) 



After simple manipulations, Y k (uS) can be rewritten as 



Y k (tu) 
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where Y^ eg (lo) is a regular function of frequency. We now 
substitute the functional Y k (ui) into Eq.(Q) for 12- The 
integral in (|]) contains a full derivative but, similarly to 
what we previously obtained in the calculation of I±, it 
does not vanish. The point is that due to a nonanaliticity, 
the logarithmic term in Eq.(|l4|) does not allow one to 
reduce the integral to just the values of Y k (u>) at the 
limits of integration. Doing the same manipulations as 
with Ji, we obtain 



I 2 = J2 ( 9 (4 + ) + ©(4) - 26(e fe )) 



(15) 



Collecting now the I\ and I2 contributions to N we 
indeed reproduce the correct result, Eq.(||). It is es- 
sential, however, that I2 is finite (e.g., for large A, 
I2 = (1 — x)Nbz )■ On the other hand, we have checked 
that if we use the perturbative expression for E, we ob- 
tain I2 in the form of an integral of a full derivative over 
frequency of a regular function Y^ ert [ui) which vanishes 
at lu = ±00. In this case, the integral is indeed equal to 
zero. 

To make this point clearer and also to discuss some 
analogy with field theory, we now compute I2 directly by 
substituting the expressions for the full Green's function 
and the full self-energy into Eq.(H). We find 



x[(uj - E+ + i5 UJ )(u - E k +iS w )Y 



(16) 



where 8 U = 5 sgnoj. In a perturbation theory, one as- 
sumes that A is small and expands I2 in powers of A: 



I 2 = £^=1 A 2 "4 n) , where l\ n> is evaluated at A = 0. 
Performing this expansion in ( |l6| ) and integrating over 
frequency, we find 



(H) 



r(n) 
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Q(e k ) - Q(e k+Q ) 
(e k -e k+Q ) 2n 



(17) 



Clearly, = for all n by symmetry: the summa- 
tion over k in ( |l7j ) goes over the full Brillouin zone while 
the function changes sign when k and k + Q are inter- 

(n) 

changed. In principle, the very fact that I 2 = to all 
orders in perturbation theory does not exclude an expo- 
nential dependence of I2 on 1/A 2 . However, we explicitly 
computed I2 numerically using Eq.(|l5|) and found that 
I2 oc A 2 at small A, i.e., I2 possesses a regular depen- 
dence on the expansion parameter. For t' = —0.2t and 
(i = -0.69t, we found I 2 ~ -0.63A 2 /t 2 . 

We now show that the reason for the discrepancy be- 
tween the explicit and perturbative calculations lies in 
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the fact that HrriA^a I 2 /A 2 =/= Irp where I^' is given 
by Eq.(p7|). Indeed, examine the form of 1% before the 
frequency integration is performed. We have 



-2i 



k 



du> 



2tt (u> - e k+Q + iS^) 2 (lj - e k + i5 u ) ' 



We see that the integrand in 1% contains a double pole 
at lo = lk+Q, i.e., the frequency integral is in fact linearly 
divergent. This divergence, however, does not show up 
in the final result as changing k to k + Q, one finds the 
same divergence, but with opposite sign. Compare now 
this formula with the exact (nonperturbative) form for I2 , 
Eq. ( |l6[ ) . We see that the inclusion of the SDW self-energy 
correction into the Green's function splits the double pole 
at lu = lk+Q into two closely located poles at u) = e k +Q 



and 



E k for ek > ek+Q, 



or uj 



for e k < £fcj 



Q- 



For most of the Brillouin zone, both poles lie in the same 
half-plane. However, for e k e k +Q < A 2 , the poles at e k +Q 
and, e.g., E~£ are located in different half-planes in which 
case the frequency integration yields a contribution which 
is absent in the case A = 0. The region in momentum 
space where this new contribution exists is very narrow, 
its area is of the order of A 2 . However, everywhere in this 
region, the denominator in 1% is also of the order of A 2 as 
Efr—lk+Q ~ A 2 . As a result, the additional contribution 
to /2/A 2 remains finite as A — > 0. To calculate this con- 
tribution, we transformed the summation over momenta 
into the integration over energies and integrated over a 
narrow region in e k+Q : \e k+Q \ < A 2 /|e fc |. We then ob- 
tained I 2 /A 2 in the form of a one-dimensional integral 
over de k : 



Urn A ^o ^2/A 2 
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where W{e k ) = ((1 + a) 2 - b 2 )(b 2 - 4a)) 1 / 2 , a = 
(c fc - 2|d)/(8|f |), b = e k /{At), e™«* = 4| M |/(1 + (1 - 
l^'l/t 2 ) 1 / 2 ), f™ = 2(4£'-/x)/(l + 2|i'|/£) < 0. Per- 
forming the one-dimensional integration numerically, we 
found for t' = -0.2< and fi = -0.69*, I 2 » -0.63A 2 /i 2 
which is exactly the same result that wc obtained using 
Eq. ( |l5[ ) . We have checked that this equivalence also holds 
for other values of i' and fi. Notice that contrary to naive 
expectations, the integration in Eq.jT^) is not confined to 
a region near the points where the Fermi surface crosses 
the magnetic Brillouin zone boundary (hot spots). These 
points correspond to e k — in the integrand in Eq.(18). 
In fact, if we substitute W(e k ) by W(0), we would obtain 
a very different I2 which for some t' even possesses the 
wrong sign. 

The anomaly which gives rise to a finite value of I2 is 
very similar to the chiral anomaly in quantum electro- 
dynamics 0]. In both cases, we have an integral which 
apparently vanishes by symmetry, but contains a hidden 
linear divergence. Performing a regularization, one gets 
rid of the divergence but simultaneously brakes the sym- 
metry. In our case, the proper regularization is achieved 



by performing all computations at a small but finite A, 
i.e., keeping the symmetry broken at all stages of the 
calculations. After regularization, the double pole gets 
split into two poles, and there appears an extra contri- 
bution from the region, where the poles are in different 
half-planes. This extra contribution turns out to be finite 
because the smallness of the phase space is fully compen- 
sated by the smallness of the energy denominator. 

Finally, let us discuss how the opening of the gap A 
influences the chemical potential /i. Combining Eqs.(||) 
and (§, we find J 2 = 2£ fe [6(e fc - M A=o) - e(e fe - fi)}. 
Clearly, the fact that I2 is finite implies that /i ^ Ha=o- 
For large A, \fi\ ra A > |/x A =o|- 

In summary, in this paper we obtained the Luttinger 
relation for the antiferromagnetically ordered state. For a 
formalism where one does not introduce an antiferromag- 
netic order parameter, we have demonstrated that the 
perturbative derivation of the Luttinger theorem breaks 
down because of a hidden anomaly which needs to be reg- 
ularized. We found similarities between this phenomenon 
and the chiral anomaly in quantum electrodynamics. 

The present paper only considers the magnetically or- 
dered state. It was recently argued by two of us || that 
I2 remains finite even in the paramagnetic phase if the 
spin-fcrmion coupling strength exceeds some threshold 
value. This argument is consistent with the idea of small 
hole pockets in heavily underdoped cuprates. A more 
detailed study of the possible violation of the Luttinger 
theorem in the paramagnetic phase is clearly called for. 
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